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We have solved Einstein's equations of general relativity for a homogeneous and isotropic met- 
ric with constant spatial curvature and found a non vanishing Weyl tensor in the presence of an 
anisotropic pressure component of the energy-momentum tensor. The time evolution of the space- 
time is guided by the usual Friedman equations and the properties of the spatial components com- 
prise a separated system of equations that can be independently solved. The physical features of 
this solution are elucidated by using the Quasi-Maxwellian equations of general relativity which 
directly connect the anisotropic pressure to the electric part of the Weyl tensor for the cosmological 
fluid. 

PACS numbers: 04.20.-q 



Oh. 
I 

O 
• I— I 

Oh' 



> 

o 

rn 

O 

m 



X 
J3 



I. INTRODUCTION 



In the literature we can find at least two formulations 
of the Einstein gravitational theory: one of them is the 
general relativity (GR) [l| and the other one is the Quasi- 
MaxwcUian (QM) equations 0. According to Lichncrow- 
icz's theorem [3|, these formulations are equivalent if ap- 
propriate initial conditions are used. It is clear that the 
theorem constrains the solutions of the Quasi-Mawellian 
equations in order that the two formulations surely give 
the same results. The main point of this restriction ap- 
pears when we want to determine the initial data in terms 
of observable quantities aiming to guarantee that the so- 
lution we shall obtain has something to do with the em- 
pirical problem we started with. It is well known that the 
observables in general relativity can only be determined 
by the geodesic deviation and that they are represented 
by the curvature tensor. This tensor has 20 independent 
components that can be expressed in terms of the Ricci 
tensor, the scalar curvature and the Weyl tensor. In the 
usual formulation of Einstein's equations only the Ricci 
tensor and the scalar curvature arc involved. In this case, 
the Weyl tensor remains completely undetermined — since 
the initial data are gij and gij^ (up to a diffeomorphism) , 
which cannot be associated to the empirical data corre- 
sponding to the Weyl tensor. In other words, it is not 
possible to look for a solution of general relativity that 
corresponds to a determined Weyl tensor specified by a 
set of empirical data. This information cannot be ex- 
pressed in terms of the initial values of the variables of 
the Einstein equations. 

All these difficulties disappear in the Quasi-Maxwellian 
formulation of gravity. In this approach, the variables 
of the theory represent directly the empirical data and 
the Einstein equations are used to relate the Weyl ten- 
sor to the energy-momentum tensor through differential 
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equations involving both. In particular, in this paper we 
present how this relation takes place when we analyze 
the standard cosmological model (SCM) closely. 

Moreover, the SCM experiences difficulties in which ex- 
otic components of matter and energy are introduced as 
an attempt to explain, for instance, the apparent acceler- 
ated expansion of the universe (dark energy), the galaxy 
rotation curves (dark matter) or the initial singularity 
in the past. Some authors claim that the main problem 
concerns the huge simplification of the geometry adopted 
and hence they suggest modifications in the spacetime 
symmetries as presented by the inhomogcneous models 
4|, in particular indicating different average processes 
5j, or modifications in the coupling between matter and 
geometry as shown by the bouncing cosmologies Q, or 
even nonconvcntional proposals as [7| , which change com- 
pletely our understanding upon the spacetime itself. 

Notwithstanding, our proposal is simpler. Instead 
of assuming modifications of geometry, we show that 
Robertson- Walker (RW) metrics admit a more general 
solution when we introduce an anisotropic pressure term 
TT^i, in the Einstein equations. Such term has not 
been considered in the establishment of the cosmologi- 
cal model by two simple reasonings: first, it apparently 
breaks the homogeneity and isotropy of the spacetime; 
second, in the case of free-shear geometries (e.g. the 
Friedman metric) there is no other tracclcss symmet- 
ric tensor phenomenologically linked to i:^u- Both ar- 
guments are not true and we shall sec that this term ex- 
emplifies exactly what wc previously pointed out about 
the problem of initial conditions in GR. 

Summarizing, in the present paper we show that if one 
considers a RW metric and an energy-momentum tensor 
of a perfect fluid with an anisotropic pressure, the QM- 
equations allow a more general solution of the Einstein 
equations which possesses, on one hand, a time evolu- 
tion given by the usual Friedmann equations and, on the 
other hand, a nonzero Weyl tensor given in terms of the 
anisotropic pressure. 
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II. FRIEDMAN SOLUTION REVISITED 



The observational data indicates that the Friedman 
geometry is more convenient to describe the homogene- 
ity and isotropy detected at large scales of the Universe. 
These symmetries suggest that the only possible fluid sat- 
isfying these properties is a perfect fluid with p (energy 
density) and p (isotropic pressure). However, as we men- 
tioned before, this is not a paradigm. In parts this is the 
scope of our paper in the sense that we will not consider 
only a perfect fluid as source for the gravitational field. 

To make a self-consistent exposition of our results, we 
present here a brief derivation of the Friedman model. 
Thus, let us start considering the infinitesimal line ele- 
ment given by 
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Alternatively, we can combine Eqs. © and ([3]) to ob- 
(1) tain the equation 



where t represents the cosmic time and a{t) is the scale 
factor. A straightforward calculation gives the following 
scalar curvature 
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where dot means time derivative and prime means deriva- 
tive w.r.t. X. The spatial curvature '■'^^i? of the hyper- 
surface defined by t = const, is 



One can sec that each value for the spatial curvature 
corresponds to a single curve in the space of solutions. 
It obviously happens because Eqs. ^ and ([7]) must be 
satisfied at the same time. Remark that solutions with 
the same sign for e are topologically equivalent. It means 
that sign{e) is enough to characterize all the solutions 
of Eq. ([6]). Therefore, these equations have only three 
relevant solutions, i.e. 
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Assuming that ^^^R has the same value everywhere, we 
set 

where e is a constant. Therefore, the scalar curvature 
becomes 



e = =^ a^x, 
e = l =^ cr = sin X, 



e = -1 



(8) 



cr = sinh^. 



We can make a coordinate transformation given by 
r = a{x) to explicitly exhibit the metric in spherically 
symmetric coordinates. The line element thus becomes 
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This equation shows that the scalar curvature of the 
spacetime depends only on time. 

The energy-momentum distribution is described by a 
perfect fluid with energy density p, isotropic pressure p 
and comoving 4- velocity = (5q , namely 

T^_,u = (P + p)Vi_yu - Pgfiy 



ds^ ^dt^ -a^{t) 



dr^ 



(9) 



Note that this metric is conformally equivalent to 
Minkowski, de Sitter and anti-de Sitter for e equal to 
0, —1 and 1, respectively. Using these coordinates, the 
time evolution of the scale factor is given by 
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In general, one assumes the existence of an equation of 
state such that p = Xp, where A is a constant. In this 
way, we can calculate the nonnull components of Ein- 
stein's equations G^^ = —T^^, (with Einstein's constant 
equal to 1), which are explicitly given by 



- = -^(l + 3A)p, 
a 6 
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where H = h/a is, called the Hubble parameter. The 
system of equations (|10p corresponds to the well-known 
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Friedman equations. A standard analysis of this model 
can be found in Q and references therein. 

In terms of the Quasi- Maxwellian formalism, the Fried- 
man equations are equivalent to 



e + - = --{i + 3X)p, 

P+(1 + A)p6i = 0, 



(11a) 
(lib) 



which are the Raychaudhuri equation and the continuity 
equation, respectively. 
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where we assumed that 



1 _/('^) 

TT 1 = 7^. 



The factor 1/a^ in the expression of irl is introduced for 
consistency reasons. Alternatively, we can combine Eqs. 
([T3| and © to obtain 



III. FRIEDMAN MODEL IN THE PRESENCE 
OF AN ANISOTROPIC PRESSURE 

In this section we basically reproduce last section 
adding to the Einstein equations an extra term on the 
right hand side corresponding to the presence of the 
anisotropic pressure tt^j^. In the formalism of fluid me- 
chanics Q, this term describes all processes involving 
viscosity and, consequently, energy dissipation. Notwith- 
standing, with the advent of the SCM in the eighties, we 
usually do not find sources for the gravitational field in- 
cluding such term anymore. At most, we can see works 
dealing with corrections to the isotropic pressure produc- 
ing a Stokesian type fluid [l3,[ll|. Concerning only the 
thermodynamics, some references analyze phase transi- 
tions produced by the gravitational field in the presence 
of the anisotropic pressure [l^ [isj . However, none of 
them take into account the anisotropic pressure on the 
right hand side of the Einstein equations. 

To do so, we propose that the most general source for 
the Friedman geometry is represented by^ 



Starting from the line clement given by Eq. ([T|), the non- 
trivial components of Einstein's equations are 
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obtain 
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^ Wc cannot add the heat flux because it breaks the isotropy 
of the RW space-time. 
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However, the solutions of Eqs. and ([Ti)) arc the same 
only if wc impose 



2k 



where k is an integration constant. This condition 
emerges from the compatibility relation of the first in- 
tegrals of both equations 
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Making the coordinate transformation given by r = cr(x), 
the line element ([T]) becomes 



dr^ 



1 
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One can see that the time evolution of the RW metric 
is still given by the usual Friedman equations (jlOap and 
(|10bp . When k = Q the Friedman model is completely 
recovered. Note that this solution is not conformally 
equivalent to Schwarzschild-dc Sitter metric. Indeed, it 
is a cosmological solution with new features. 



IV. QUASI-MAXWELLIAN EQUATIONS AND 
RW METRICS 



In this section, we use the Quasi-Maxwellian equations 
of gravity to derive the solution expressed by Eq. (|16p . 
It is convenient to assume an observer field as = 
(5q . The energy momentum distribution is given a perfect 
fluid with anisotropic pressure ■n^y. Using the metric 
represented by Eq. ([T|). these assumptions lead to the 
following non vanishing Quasi-MaxwcUian equations (see 
details in Appendix): 
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p + {p + p)e = 0, 
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(17a) 
(17b) 
(17c) 
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As we said before, Eqs. (|17ap and (jl7b[) correspond to the 
Friedman equations. Eq. (|17cp yields the electric part of 
the Wcyl tensor in terms of the anisotropic pressure. The 
equation (jl7d[) represents the compatibility condition of 
the Einstein equations with the (constant) spatial curva- 
ture given by Eq. ([15]). Remark that this equation does 
not imply that the Weyl tensor is identically zero. 

From a straightforward calculation, the electric part of 
Weyl tensor is read as 



1 







(18) 



The first integral coming from the (constant) spatial 
curvature equation is expressed by 
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The common solutions are obtained by setting Ci = 
and C2 = —2k. 

Finally, we make the coordinate transformation r 
a{x) to rewrite the line element ([1]) as given by Eq. ([1 
Note that this solution has a non vanishing electric part 
of the Wcyl tensor given by 



-10 
1/2 
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and the spatial curvature ^'^^i? remains constant. Remark 
that Eq. (P5|) is very similar to the Newtonian tidal forces 
times a time dependent function 1/a^ . The consequences 
of this fact will be discussed in next sections. 



TRAJECTORIES OF TEST PARTICLES 



where 
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Substituting this expression into Eq. (jl7d[) . we have to 
solve an extra constraint, which is given by 
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Instead of integrating completely the geodesic equa- 
tions of the metric ([T6l) , which seems unnecessary at this 
moment, a qualitative analysis of the particle trajectories 
in this geometry can be done. In fact, we only analyze the 
geodesic motion of test particles (timelike and nulllike) 
moving along the equatorial plane {9 = 7r/2 and ^ = 0). 
This simplification allows to study the behavior of the 
effective potential to which these particles are subjected 
and comparison to other cases can be easily done. 

In this way, the simplified geodesic equations reduce to 



where h[t) is an arbitrary function. Multiplying both 
sides of Eq. (pil)) by a^cr^, the following relation appears 



2 II /2 I 



-ia^{t)h{t). 



(21) 



aa 



A 



(26) 
(27) 



The equation for the time coordinate is trivially satisfied 
by fixing h = —kja?^ where fc is a constant. To solve the 
equation for the spatial coordinate x, we use Eq. the 
spatial curvature in terms of cr, and rewrite Eq. (j2ip as 
follows 



2-i'. 
a 



2-6' 



f = 0, 



(28) 



where we denoted X' = dX/dr, t is the affine parameter 
along the curve (which is the proper time for the timelike 
geodesies) and b is equal to 1 for timelike geodesies or 
for nulllike ones. We also denoted 



0, 



(22) 



Multiplying both sides of this equation by 2(Ta' , we get 
a first integral of this form: 



a 



(23) 



A{r) 



2k 

r 



Remark that we do not fix the value of the spatial cur- 
vature in the analysis below: e remains an arbitrary con- 
stant. 

First, we solve Eq. (|28p and see that the angular mo- 
mentum is a conserved quantity 
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where / is an integration constant. Then, substituting 
Eqs. (HH) and (EH) into Eq. ^ we get 

t" + ^(t'2_^) = 0, (30) 
a 

which can be integrated, yielding 

it"-b) = -^, (31) 

where E is another integration constant, in general asso- 
ciated to the total energy of the test particle^. Substi- 
tuting this equation in Eq. (^7)) results 

a\'' = - (^E + A. (32) 

This equation can be seen as the energy conservation 
equation of a particle moving in a one-dimensional effec- 
tive potential. Note that this equation has a constant 
term E (like a mechanical energy), a kinetic- like term 
a^r'^ and the remaining ones correspond to the effective 
potential that we denote by V{r). 

To compare Eq. ([5^ with the effective potential ob- 
tained from the Schwarzschild metric we suppose 
the spatial curvature equal to zero (e = 0). Therefore, 
Eq. ^ becomes 



where we seek for the planetary orbits u = u{(j)). How- 
ever, the time-dependence of Eq. ([5^ makes the problem 
more complicated. Therefore, we consider an interval of 
cosmological time (<07^i) in which the scale factor re- 
mains almost constant, i.e., a(t) = oq for t £ (io,ii)- 
In other words, we assume that the cosmological time 
changes very slowly in comparison to the period of the 
planetary orbits. We thus obtain 

f (—X = -lE + Pu^)ll-2ku). (34) 
\d(t)J 

Differentiating this equation with respect to 0, yields 

-—iT+u = 6ku + (35) 
aq) L 

This is exactly the differential equation given by the 
Schwarzschild metric in response to the Kepler problem 
if we identify k to the effective mass of the gravitational 
source and E is identified to the total mass of the test 
particle. In particular, considering i? = 1 we can cor- 
rectly reproduce the perihelion shift and for i? = we 
obtain the light deflection predicted by general relativity 
via Schwarzschild metric. Note that 6 and E practically 
coincide in this regime. In this way, we suggest that dark 
matter effects can be interpreted as a convenient choice 
of the Weyl tensor as shown by Eq. pS]) or even that 
dark matter is only a problem related to the bad choice 
of initial conditions in GR. 



-I ,o „ QikE I '2ikt 
a\'^ ^ -E + — - — + — . (33) 

The effective potential V{r) can be written down as 

2kE P 2kl^ 
Vir) = + - - — . 

Note that this potential presents the same terms as the 
Schwarzschild geometry in the case of a single particle 
moving along a geodesic in this metric. In particular, 
both potentials have the same qualitative behavior if we 
choose k> Q and E > Q and, therefore, they produce the 
same particle trajectories. 

For the sake of completeness and comparison, we start 
a more detailed, but still qualitative, analysis of the Ke- 
pler problem in this solution rewriting Eq. (|32p in terms 
of the variable 




^ Opposite to the Schwarzschild case, the first integral of the time 
component of the geodesic equation is not necessarily positive 
definite. 



VI. GEODESIC DEVIATION OF THE 
COSMOLOGICAL FLUID 

Among its attributes, the Quasi-Maxwellian represen- 
tation of gravity has the quality to put together the for- 
malism of the electromagnetic interactions with a for- 
mal approach of GR. Nonetheless some fundamental dis- 
tinctions must be stressed. The empirical determination 
of an electromagnetic field is made through the Lorentz 
force and a test particle in order to identify the pres- 
ence of an electromagnetic field. The electromagnetic 
tensor, obtained from the integration of Maxwell equa- 
tions, does not distinguish the contribution of the local 
charges and currents distribution from boundary condi- 
tions. In GR the empirical identification of a gravita- 
tional field cannot be done using a single test particle 
since the Christoffel symbols can be set equal to zero 
by coordinate transformations. In order to empirically 
determine the properties of the gravitational field it is 
necessary to look for the geodetic deviation expressed in 
terms of the curvature tensor. This tensor explicitly sep- 
arate the contribution coming form the local distribution 
of the energy-momentum tensor, algebraically associated 
to the traces of the Riemann tensor, from the global con- 
tribution of boundary conditions represented here by the 
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Weyl tensor. Indeed, the measurements of the gravita- 
tional field effects can be only done through the geodesic 
deviation equation which determines the rate of the rel- 
ative acceleration between infinitesimally near geodesies. 
This equation is then given by 



(36) 



where is the deviation vector and V^^ is the vector 
field tangent to the geodesic congruence. 

The distortion produced by the Weyl tensor upon a 
given congruence of curves can only be detected by this 
equation substituting the Ricmann tensor by its decom- 
position into irreducible parts: the Ricci tensor, scalar 
curvature and the Weyl tensor. 

Since we are dealing with the cosmological fiuid we set 
= (5q , evaluate the right hand side of Eq. ([55)1 and 

get 



Now comes a remarkable result: according to Eq. (jl7cp 
the term inside big brackets is identically zero. There- 
fore, the cosmological fluid does not measure any distor- 
tion caused by the presence of the anisotropic pressure. 
In other words, the distortion caused by the anisotropic 
pressure and the electric part of the Weyl tensor are com- 
pensated in such a way that the cosmological observer 
does not attribute any eventual modification of the space- 
time to these quantities, enabling one to set iJ^j, and tt^^ 
equal to zero by hand. However, this is not allowed if 
we want to understand correctly the gravitational field 
effects in the Universe using the empirical data as initial 
conditions. 



VII. CONCLUDING REMARKS 

We could see that general relativity does not contain in 
its dynamics all the information necessary to determine 
the curvature tensor from empirical data. It means that 
the initial condition problem in GR (also called Cauchy 
problem) should be revisited, since a more realistic de- 
scription of the Universe must be privileged (an open 
system rather than a close totality). 

In particular, in this paper, we showed how the stan- 
dard cosmological model sets the Weyl tensor equal to 
zero ah initio and that it is not consequence of the 
Friedmann equations. Therefore, we have developed a 
cosmological model with constant spatial curvature and 
nonzero Weyl tensor without wasting homogeneity and 
isotropy. 



VIII. APPENDIX: QUASI-MAXWELLIAN 
EQUATIONS 

Bianchi's identities together with the Einstein equa- 
tions of gravity yield the Quasi-Maxwellian equations if 
we rewrite the Riemann tensor in terms of its traces and 
the Weyl tensor. Assuming Einstein's constant equal to 
1, these equations are 



W 



1 



"^2 6 
The Weyl tensor is given by 



1 



(37) 



where the auxiliary tensors are 

2Maj3fiv — Ra^igpu + RpyQaii — RavQp^i — Rp^iQav 

and 

9ap^iv = gan9pv — 9au9i3fi- 

The reason of this nomenclature is due to several analo- 
gies between the Quasi-Maxwellian and the Maxwell 
equations. However, this similarity does not go far be- 
cause QM-equations are in fact highly non-linear in com- 
parison to Maxwell's theory, leading to situations that 
never happen in the last case. 

The projection of the QM-equations with respect to 
the vector field V" and its orthogonal hypersurface al- 
lows us to do calculations with. At this point, it is very 
useful to replace the Weyl tensor by its electric and 
magnetic H^ia parts, respectively: 

where *Wafipu is the dual of Weyl tensor constructed 
using the skew-symmetric Levi-Civita tensor. 

The covariant derivative of can be decomposed into 
its irreducible parts: 



where 9 = V^-fj. is the expansion coefficient, = V^-^V^ 
is the acceleration. 



1 . 



is the shear tensor and 
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is the vorticity. Therefore, four independent projections 
lead to the following linearly independent equations 

(38) 
(39) 

(40) 

W'^^v'^^'^v^VrE^^e^p - ^Hpt^.^hi'T^^^^^Pv^ = 

J/i^^(g^;^ - qt^a,, - TT'^^a^,) - ^(p + pW'^ + ^g(^a^) + 

(41) 

It is easy to see the similitude to Maxwell's theory. 
The conservation law of the energy-momentum tensor 
T^'^;^ = gives 

p+ip + p)e + + q"-c.- tt'^V^, = 0, (42) 



(p + p)aa - Ptj,h^a + qfih^a + 0q 

(43) 

+q''LJai. + TTa";^ + Tr'"'<J^^Va = 0. 

The integrability condition can be translated into 
equations for the kinematical quantities. Thereby, these 
evolution equations are 



9 + - + 2{a' + uj-')-a",o. = -^{p + ip), (44) 

-'■2ha'^hi3''{af,.^ + au■^J,) + ^6aap + crap(T^/3+ (45) 

(46) 

together with the constraint equations 



^e,p/i^-(^"7+^"7);a^^-«"K.+c^A.) = -qx, (47) 



+2w"aa = 0, (48) 



iJrA = -^/l(r'/lA)"?7£'^^''K(cra/3+a;„/3);^+a(^tJA). (49) 

This set of equations corresponds to the necessary 
equations to propagate GR on a given hypersurface to 
the whole spacetime. 
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